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ABSTRACT 

We use a temperature map of the cosmic microwave background (CMB) obtained using the South Pole Tele- 
scope at 150 GHz to construct a map of the gravitational convergence to z ~ 1100, revealing the fluctuations in 
the projected mass density. This map shows individual features that are significant at the ~ 4a level, providing 
the first image of CMB lensing convergence. We cross-correlate this map with Herschel /SPIRE maps cover- 
ing 90 square degrees at wavelengths of 500, 350, and 250 |Ltm. We show that these submillimeter-wavelength 
(submm) maps are strongly correlated with the lensing convergence map, with detection significances in each 
of the three submm bands ranging from 6.7 to 8.8 a. We fit the measurement of the cross power spectrum 
assuming a simple constant bias model and infer bias factors of b = 1.3 — 1.8, with a statistical uncertainty of 
15%, depending on the assumed model for the redshift distribution of the dusty galaxies that are contributing 
to the Herschel /SPIRE maps. 

Subject headings: galaxies: structure — cosmic background radiation 
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1. INTRODUCTION 

Gravitational lensing of the cosmic microwave background 
(CMB) is emerging as a powerful cosmological tool. The spa- 
tial variation of the statistical properties of the CMB that is 
induced by gravitational lensing was first detec ted in cross- 
correlation with radio- sel ected galaxy catalogs ( Smith et al. 
2007 [ Hirata et al. 2008| ), and subsequently detected inter- 
nally in CMB maps by the Atacama Cosmology Telescope 
collaboration (Das et al .|2011|) and the South Po le Telescope 
(SPT) collaboration (|van Engele n et al.|2012| . In Bleem et al. 
(2012), reconstructions of the mass distribution were found 
to correlate strongly with galaxy cat alogs selected in both 
the optical and infrared bands, while Sherwi n et al. (2012) 
showed that CMB lensing was well-correlated with quasars. 

Using the CMB as the background source to study grav- 
itational lensing by intervening structure offers several ad- 
vantages over using distant galaxies: the source redshift is 
the same for all lines of sight, is extremely well-known, and 
has the highest redshift observable with electromagnetic ra- 
diation. The statistical properties of the source are well- 
characterized, and CMB maps cover areas ranging from a few 
hundred square degrees to the full sky. However, the single 
redshift for the CMB does not allow any information about 
the redshift distribution of the mass along the line of sight, 
and noise levels in CMB lensing convergence maps at cur- 
rent sensitivities are substantially higher than noise levels in 
cosmic shear measurements. 

As CMB lensing is an integral along the entire line of sight, 
the strongest cross-correlations will be with sources that have 
a similarly broad extent in redshift sp ace. As demons trated 
below, and as theoretically predicted (Song et al. 2003]), the 
cosmic infrared background (CIB) fluctuations provide an ex- 
cellent match. The CIB at submm wavelengths is believed to 
have a substantial contribution from sources from redshifts 
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In this paper, we cross-correlate a map of the gravitational 
lensing convergence (proportional to the surface density) de- 
rived from SPT temperature data at 150 GHz with maps of the 
submillimeter-wavelength (submm) sky at 500, 350, and 250 
[im obtained with Herschel /SPIRE. By using maps rather 
than catalogs, as was done in previous CMB lensing cross- 
correlations, we are studying emission from sources that are 
individually unresolved. The SPT and Herschel datasets are 
described in sections 2 and 3, and the results of the cross- 
correlation are presented in section 4. A comparison with a 
simple theoretical model is presented in section 5, and we 
conclude with a discussion of the results. 

2. CMB MAP AND CORRESPONDING MASS MAP 

The SPT has been used to image 2500 square degrees to a 
depth of <18 /iK-arcmin at 150 GHz, and two ~ 100-square- 
degree fields (each subtending lh in Right Ascension and 10 
degrees in Declination) within this area to a depth of ~13 
/iK-arcmin. For this work, we use observations centered on 
one of those deeper fields, centered at (RA,DEC)=(23h30m, 
-55d00m), using data from both the 2008 and 2010 observing 
seasons; the recen t CMB power spectrum measurements of 
|Story et al.| ( [20T2| ) used only the data from 2008 for this field . 

A CMB map is generated as outlined in Story et al. ( 2012| ). 
In addition, to avoid apodization effects at the edges of the 
field when constructing the lensing map, data from surround- 
ing fields are used to make a single larger CMB map 17.1 



degrees on a side. This map extends well beyond the re- 
gion covered by Herschel data. The input CMB map is shown 
in the left panel of Figure [T| Adjacent fields are combined 
using inverse- variance weights in overlapping regions, and 
there is no evidence for any discontinuities at the boundaries. 
Point sources and massive galaxy clusters are removed usi ng 
a Wiener-interpolation algorithm ( van Engelen et al.|2012 ). 

Simulated CMB maps are obtained by coadding simulated 
signal and noise realizations for each individual SPT field. 
The simulated maps are made with known input gravitational 
potentials, and simulated signal maps are generated using 
timestream-based simulations, as in |Story et al.|p012] ). Noise 
realizations are obtained directly from the observations, by 
taking randomized combinations of the data which remove all 
sky signal, as detailed in van En gelen et al.| ( |2012] ). 

The analysis procedure is applied to both the real and sim- 
ulated SPT maps. Gra vitational converge n ce ma ps are gen- 
erated as outlined in |van Engelen et al.] (|2012|), using th e 
quadratic estimator method ( |Hu|2001[|Hu & Okamoto(2 002 ). 
This method entails constructing a gradient-filtered map and 
an inverse- variance weighted map (i.e., two different filterings 
of the same CMB field), multiplying them together and tak- 
ing a divergence. The resulting product can be shown to be an 
estimator for the map of the gravitational potential. The effec- 
tive transfer function due to the SPT filtering was constructed 
by cross-correlating the derived lensing potential of the sim- 
ulated maps with the lensing potential maps used to generate 
those simulations. 

Foreground contamin ation of the lensing converg ence maps 
is expected to be small: van Engelen lk al.| ( |2012] ) found that 
residual contamination of the lensing convergence map from 
point sources and galaxy clusters is expected to be at the level 
of a few %. The sign of this effect is expected to be neg- 
ative, such that foreground contamination acts to reduce the 
observed cross-correlation. 

The resulting lensing convergence map is shown as con- 
tours in Figure [2] Features can be seen with significances 
exceeding 4a. 

3. Herschel /SPIRE MAPS 

Submillimeter maps at 500, 350, and 250 ja m are created 
using observations with the SPIRE instrument ( Griffin et al. 



2003| aboard the Herschel Space Observatory (Pilbrat t"et al. 
2010) obtained under an OT1 program (PLCarlstrom). Obser- 
vations were made in SPIRE fast-scan mode (60 arcsec s _1 ) 
and consisted of two sets of orthogonal scans covering 
~ 90 deg 2 . The observing strategy was chosen to optimize 
sensitivity to large-scale signal and provide redundancy for 
measuring the auto-frequency power spectrum of background 
fluctuations. 

Maps are made with SMAP, an iterative mapmaker designed 
to optimally separate large-scale noise from signal; the map- 
makin g algorithm is d escribed in detail i n |Levenson et al.| 
P010| ) and updated in |Viero et"aL] ( |2012| ). To estimate the 
transfer function we use the same map-making process on 
mock SPIRE data. For both real and mock data we make maps 
with 10 iterations; we have checked that the maps are ade- 
quately converged at this point. Additionally, time-ordered 
data (TODs) are divided into two halves and unique "jack- 
knife" map-pairs are made. To avoid having to reproject 
or regrid the Herschel /SPIRE maps, we make them using 
the Lambert azimuthal equal-area projection (also known as 
zenithal equal area, ZEA), with astrometry identical to that of 
the SPT map, and with 30" pixels. 
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FIG. 1 . — SPT 150 GHz (2mm) temperature map (left) and Herschel /SPIRE maps (right) used for this analysis. For display purposes only the inner *• 
the SPT temperature map that was used to construct the lensing map is shown. In the right panel, (red, green, blue) correspond to (500, 350, 250) |om. 

The maps have rms instrument noise levels (per 30" pixel) 
of 14, 10, and 7mJy, while the instrument effective point- 
spread functions are 36.6, 25.2, and 18.1 " full widths at half 
maxima (FWHM) at 500, 350, and 250 (xm, respectively, The 
30" pixelization of the maps reduce the resolution substan- 
tially on small scales, but pixelization and instrument noise 
effects are not important on the scales of int erest for this study. 
The last step, following |Viero et aL] ( |2012| ) is to convert the 
maps from native units of Jy beam 1 to Jy sr _1 , which is 
done by dividing them by the effective beam areas, 3.688, 
1.730, and 1.053 x 10 -8 steradians. Color corrections from a 
flat-spectrum point-source calibration have a negligible effect. 
The absolute calibration is accurate to 7%, an uncertainty that 
is small compared to our statistical precision. 



< 60% of 



sume a Gaussian error distribution. Cross powers are reported 
in Table [T] and are shown in Figure [3] 

The signal-to-noise ratio in the cross correlation is sub- 
stantial: at 500, 350, and 250 |Ltm the model with no cross- 
correlation is strongly disfavored relative to the best-fit lens- 
ing amplitude, with \ 2 differences of 79, 69, and 45, respec- 
tively. Lensing is detected in every power spectrum bin. 

5. THEORETICAL MODEL 

As a cross-check on the shape and amplitude of thes e spec 
tra, we adop t the simple constant bias model used in Bleem 
et al. (2012), using the non-linear power spectrum at each red 
shift: 



4. RESULTS AND ANALYSIS 

In Figure [2] we present convergence and submillimeter- 
wavelengths maps filtered to emphasize modes in the lensing 
map that have significant (>0.5) signal-to-noise, allowing a 
by-eye comparison of the structure. Modes with L < 100 
(scales larger than 2°) have been filtered to remove scales 
where the timestream filtering of the submm-wave maps be- 
comes substantial. The SPT temperature map has spatially 
anisotropic noise (Schaffer et al.|[20TT] ), whi ch ultimately 
leads to an isotropic noise in the lensing map ( van En gelen 
et al.||2012| ). This leads to a tendency for modes to be bet- 
ter measured when they have more horizontal structure than 
vertical structure. 

Due to the imperfect redshift overlap, the lensing map can 
have features that are not in the submm maps; in particu- 
lar, high-redshift structure (z > 3) will appear relatively 
stronger in the lensing map, while structure below z ~ 0.5 
will be strongly suppressed in the lensing map as compared 
to the submm map. The submm maps are extremely well- 
correlated with each other, while the lensing map has several 
features that are not well-matched in any of the submm maps. 
Nonetheless, there are many features in common between the 
maps. 

To compar e these maps quantitatively, we use cross-power 
spectra, as in |Bleem et al.| ( [20T2| ). Uncertainties are obtained 
by cross-correlating each submm map with lensing mass maps 
obtained from simulated SPT maps. We use the rms in cross- 
power simulated amplitudes as the rms uncertainty and as- 



C£=bj dz^^W*(x)W l (x)PDM 



(1) 



where W K (\) gives the redshift weighting of the mass map 
and W l (x) is propor tional to the line of sight distribution o f 
the intensity dl/dx ( [Bleem et aL||2012[ |Song et al |I2003| ). 
The non-linear power spectrum of the dark matter, Pbm, is 
calculated using CAMB and Halofit, assuming the best-fit 
WMAP9+SPT cosmological parameters for a flat ACDM cos- 



table l 

CMB CONVERGENCE-SPIRE CROSS POWER SPECTRUM 



L 


£<500 


£<350 


£>250 




(mJy/sr) 


(mJy/sr) 


(mJy/sr) 


150 


38 ± 18 


101 ± 29 


134 ± 58 


250 


49 ±21 


65 ±30 


74 ±45 


350 


10 ± 11 


18 ± 15 


26 ±21 


450 


19.0 ±6.2 


38 ± 10 


52 ± 13 


550 


8.5 ±6.4 


15 ± 10 


8± 12 


650 


13.7 ±4.6 


16.0 ±7.8 


16.9 ±8.8 


750 


13.3 ±4.2 


17.5 ±5.6 


14.3 ±8.5 


850 


4.9 ±3.2 


3.5 ±5.4 


13.5 ±7.0 


950 


6.1 ± 1.9 


9.5 ±3.5 


7.7 ±5.0 


1050 


6.9 ±1.9 


8.1 ±3.2 


2.1 ±4.6 


1150 


1.2 ±1.5 


2.9 ±2.4 


0.7 ±3.6 


1250 


5.7 ±2.2 


9.2 ±3.1 


11.6 ±4.2 


1350 


2.3 ±1.5 


5.5 ±2.6 


6.1 ±3.5 


1450 


4.2 ±1.8 


5.1 ±3.1 


3.8 ±4.1 


1550 


2.2 ±1.6 


2.4 ±2.9 


3.1 ±3.3 
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Fig. 2. — Map of the CMB lensing convergence measured with SPT data (contours in all panels) and overlaid on maps of the 500, 350, 250 jam Herschel /SPIRE 
data (top right, bottom left, bottom right, respectively). All maps have been filtered to only show scales in the lensing map that are expected to have typical signal 
to noise of at least 0.5, which suppresses all features on scales smaller than ~ 0.5°. All maps have been masked by the SPIRE coverage. Lensing contours are 
spaced by 1 a of noise. Red (blue) indicates regions of increased (decreased) mass or flux. 



mology ( Stor y et al.|2012| ). 

The redshift distribution of contributions to the submm 
background has been extensively studied in recent years, and 
there exist substantial disagreements between authors. We 



adopt t wo r ecent determ i nation s, presented in Bethermin et al. 



pOTT]) and |Viero et aLj ( |2012| ) that roughly bracket expecta 
tions, to predict the cross-correlation signal. To derive this 
signal, we assume that the submm light traces the non-linear 
dark matter density field at every redshift, with a single ampli- 
tude, the bias b, that we fit to the data. The cross-correlation 
will be most sensitive to redshifts z ~ 0.5 — 3, with lower z 
a poor match to CMB lensing, and higher z not having sub- 
stantial submm emission. As seen in the insets of Figure |3j 
the 500 |Ltm emission is expected to have broader overlap with 
the CMB lensing kernel, and should therefore show a stronger 
correlation. 



Fits are performed using points between L = 100 and 
L = 1600, as done in previous SPT lensing studies. The best- 
fit bias parameters for each observing wavelength and redshift 
distribution choice are shown in Table [2j with best-fit bias pa- 
ramete rs depending on which re dshift distribution is assumed. 
For the |Bethermin et al.| ([2QTT| ) model we find 6 ~ 1.8 ± 0.3 
while the Vier o et al.| ( [2Q12| ) model for the CIB intensity gives 
b ^ 1.3 ±0.2. The uncertainties reflect statistical uncertain- 
ties only, and the large difference between the two models in- 
dicates that systematic uncertainties are substantial. The dif- 
ference in bias factors is largely due to the different integrated 
mean intensities in the two models; for example, at 500 |j,m 
the two models predict mean intensities that differ by a factor 
of 1.5, while the derived bias factors differ by a factor of 1.4. 
This difference in the mean intensity is larger t han the ~ 25% 
uncertainty in the FIRAS measurements (Fixs en et al.|1 998 ); 
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TABLE 2 
Fits to constant bias model 



Wavelength 



Bias (V12) 



Bias (Bll) 



500 jam 1.29 ± 0.16 (12.6) 1.80 ± 0.22 (12.7) 
350 urn 1.35 ± 0.17 (9.7) 1.82 ± 0.24 (9.9) 
250 yim 1.34 ± 0.23 (11.8) 1.56 ± 0.27 (12.0) 
Cross-spectrum best-fit amplitudes to constant bias model for Viero et al 
(2012) (V12) and Bethermin et al. (2011) (Bll) redshift distributions, x 2 
of fit shown in parentheses. Quoted uncertainties only include statistical 
uncertainty. 



the mean intensity in the Viero et aT] ( [2012 ) model is more 
than 2<T higher than that measured by FIRAS at 500 |Ltm. 

This simple model provides a very good fit, with y 2 = 12.6 
or 12.7 for 14 degrees of freedom at 500 |Ltm, depending on 
the assumed redshift distribution of the submm background. 
Despite the qualitative difference in the two redshift distri- 
butions apparent in the insets of Figure |3j good fits are ob- 
tained for both models, although a different normalization is 
preferred by each. This arises because most of the power is 
coming from the non-linear regime, where a powe r-law is a 
remarkably good fit to the clustering power spectra ( Addison 



et al. 2012). As the cross-spectrum is a superposition of sim 
ilar power-laws from different epochs, the detailed redshift 
distribution does not affect the shape of the cross-spectrum. 

The bias factors at infrared and submm wavelengths have 
been measured using both source catalogs and auto power 
spectra of the diffuse backgroun ds, as reviewed re cently for 
the cosmic infrared background in Penin et al. (2012). The in- 
ferred bias values depend on the assumed redshift distribution 
and intensity of the background, and the bias value that we 
measure is the clustering amplitude relative to the non-linear 
matter power spectrum, rather than either the linear matter 
power spectrum or a halo model, so a direct comparison is 



difficult. Using BLAST data at 500, 35 0, and 250 \m\£ 
et al. 2009| ) and intensity estimates from Lagache et al. |([2' 



typical ^EIas fac t ors o f 2.2 ±0.2 were found (Penin et aL 



Viero 
•004), 



2012). 



Amblard et al.| (|2011|) find slightly higher bias values using a 



halo model and fitting internally for the intensity as a function 
of redshift. The bias values found here are somewhat lower, 
but could be explained by differences in the assumed mean 
intensities and their redshift distributions. 

Some studies of dusty source s at high redshift have led to 
substantially higher bias factors: Brodwh Tet al.| (2008 ) found 
that z ~ 2 dusty, obscured galaxie s selected in tr i e optic al/IR 
had bias factors b ~ 3 — 5, while Hicko xet al . | ( |20 1 2"| ) used 
sources selected at 870 |Ltm to estimate 6^3. 

For comparison with lower redshift galaxy samples, recent 
results from SDSS-III find bias factors of ~ 2 for the massive 
galaxies (halo masses ~5 x 1O 13 /i _1 M ) bei ng targeted for 
baryon acoustic oscillation studies at z ~ 0.3 (Parejko et al. 
|2013| ), while bias estimates based on the SDSS main galaxy 
sample ( [McBride et al.|2011| ) find b = 1 - 1.2 for typical lu- 
minosity (L*). This suggests that the typical contributors to 
the submm background could be the higher redshift precur- 
sors to (or at least have the same mean bias as) galaxies that 
are intermediate in mass between these two samples. 

In work that is closely related to the current work, Hilde- 
brandt etaL] ( |2013| ) cross-correlated gravitational lensing of 



Lyman-break galaxies with a catalog of sources detected at 
250 |Ltm, and inferred typical masses of 1.5 x 1O 13 M for 
these galaxies. 



6. DISCUSSION AND CONCLUSIONS 

We have shown that large-scale structure traced by submm 
sources is well-correlated with a CMB lensing convergence 
map. The cross-correlation is highly significant at 500, 350, 
and 250 \im, corresponding to detection significances of 8.9, 
8.3, and 6.7 a, respectively. 

The cross-correlation between the lensing convergence map 
and each submm map is well fit by a simple constant bias 
model, with bias factors of b = 1.3 — 1.8, depending on the 
assumed redshift distribution for the submm intensity. The 
lower bias factors are found for an assumed intensity distribu- 
tion with more flux coming from higher redshifts. 

There are several ways to extend the utility of the lens- 
ing convergence- SPIRE cross-power spectra presented here. 
For example, combining them with the cross-power and auto- 
power spectra among the three SPIRE bands will probe the 
redshift distribution of the contributing sources and the corre- 
spondence between submm flux and the underlying dark mat- 
ter distribution. 

This technique is highly complementary to studies of the 
auto- and cross-correlations of submm background maps. 
While convergence maps have more noise (at current CMB 
map noise levels), concerns about Galactic cirrus or separat- 
ing shot noise are greatly reduced, making cross correlation 
with CMB lensing an extremely robust probe of clustering 
with a promising future. 

With the release of Planck maps covering a broad range 
of CIB wavelengths with well-matched angular resolution, it 
will be possible to perform a similar analysis over the entire 
2500 square degree SPT survey area, while the coming Dark 
Energy Survey (DES) will also have nearly complete overlap 
with this area. DES will have both galaxy catalogs and cos- 
mic shear maps with some resolution in the line of sight direc- 
tion. In combination with the SPT CMB lensing convergence 
this will enable 3D mass maps of the universe extending to 
z~ 1100. 
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